Thermal Entanglement in Ferrimagnetic Chains 
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A formula to evaluate the entanglement in an one-dimensional ferrimagnetic system is derived. 
Based on the formula, we find that the thermal entanglement in a small size spin-1/2 and spin- 
s ferrimagnetic chain is rather robust against temperature, and the threshold temperature may 
be arbitrarily high when s is sufficiently large. This intriguing result answers unambiguously a 
fundamental question: "can entanglement and quantum behavior in physical systems survive at 
arbitrary high temperatures?" 
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A physical system may exhibit entanglement at a finite 
temperature \^, 2.i ^j As M ■ The thermal entanglement 
always vanishes above a threshold temperature for sys- 
tems with finite Hilbert space dimension q. Recently, 
Ferreira et al. raised a fundamental question: can entan- 
glement and quantum behavior in physical systems sur- 
vive at arbitrary high temperatures? Q They found that 
the entanglement between a cavity mode and a movable 
mirror does occur for any finite temperature. This result 
sheds a new light on the question and help to understand 
macroscopic properties of solids. 

In this report, we derive a formula to evaluate the en- 
tanglement in an one-dimensional ferrimagnetic system. 
Intriguingly, we find that the entanglement is rather ro- 
bust against temperature in this kind of system, with the 
Hamiltonian 



H — J 2_^{^n • S„+i + S„+i • S„+2), 



(1) 



where s„ and S„ are spin-1/2 and spin-s operators, re- 
spectively. The antiferromagnetic exchange interactions 
exist only between nearest neighbors, and they are of the 
same strength which is set to unity (J = 1). Physically, 
the system contains two kinds of spins, spin i and s, al- 
ternating on a ring (or a chain with the periodic boundary 
condition) . 

Let us now study the entanglement of states of the 
system at thermal equilibrium described by the density 
operator p{T) = exp{—(3H)/Z, where f3 = l/ksT, ks 
is the Boltzmann's constant, which is assume to be 1, 
and Z = Tr{exp(— /3i/)} is the partition function. The 
entanglement in the thermal state is referred to as the 
thermal entanglement. 

To study quantum entanglement in the ferrimagnetic 
system, we need a good entanglement measure. One pos- 
sible way is to use the negativity pj based on the partial 
transpose method 9]. In the cases of two spin halves 
and the (1/2,1) mixed spins, a positive partial trans- 
pose (PPT) (or the non-zero negativity) is necessary and 
sufficient for separability (entanglement). Although the 
present ferrimagnetic system is a kind of (1/2, s) system. 



fortunately, it was shown that due to the SU(2) symme- 
try in the model Hamiltonian (^), the non-zero negativity 
is still a necessary and sufficient condition for entangle- 
ment between a spin half and spin s |lO|. This result 
allows us to exactly investigate entanglement features of 
our mixed spin systems. 

The negativity of a state p is defined as 



j^{p) = Yl I'"* 



(2) 



where pi is the negative eigenvalue of p^^ , and T2 denotes 
the partial transpose with respect to the second system. 
The negativity N is related to the trace norm of p'^^ via 



A/'(P) = 



IIP^^IIi-1 



(3) 



where the trace norm of p"^'^ is equal to the sum of the 
absolute values of the eigenvalues of p^^ . 

Obviously, our system has the SU(2) symmetry, and 
any two-spin reduced density matrix from the thermal 
state is also SU(2)-invariant. Now, we consider the en- 
tanglement between the spin half and spin s, and derive 
the corresponding expression of negativity by the par- 
tial time reversal method [ll| , which is equivalent to the 
partial transpose method up to a local unitary operator. 

The density matrix of an SU(2)-invariant state for the 
spin half and spin s can be written in the form 



^ 2s 



-1/2 



l-F. 

2s + 2" 



s+1/2 



(4) 



with F G [0, 1]. One may immediately check that the pa- 
rameter F is identical to the expectation value of the pro- 



jector Ps-1/2 on the densitry matrix, i.e., F 



Noting the fact that Ps-1/2- 
density operator as 



-1/2 



s+1/2 



1, we rewrite the 



{2sF + F ~s)F, 



-1/2 



2s(s + 1) 



1 -F 

2s + 2' 



(5) 



and obtain the relations 



1 



s-l/2 



s+1/2 



2.S+1 
1 

'2s + l 



(s - 2si • S2) 

(s + l + 2si-S2). 



(6) 



The partial time reversal operator T2 changes the sign of 
S2 and Si • S2 



''"2(82) = ^ S2 
r2(si • S2) =-Si -82, 



Therefore, we get 

T2(Ps-l/2) 



2s 



2s +1 



s-l/2- 



(7) 



(8) 



We see that the partial time reversal also changes the sign 
of projector, but with an additional additive constant. 

Using Eqs. © and (|SJ|, the partial time reversed den- 
sity matrix is obtained as 



P 



{2sF + F-s)P,_i/2 2sF + F + l 



2s(s + l) 



2(2s+l)(s + l)' 



(9) 



Any projector has two eigenvalues and 1. Thus, from 
the above equation, we deduce that only the following 
eigenvalue of p'^^ is possibly negative 



A 



1 



2s + 1 
1 



F 

2^ 



2s 



i + i/si-Sa; 
1 V 2 s ^ ' 



(10) 



where we have used the first equality in Eq. ©. By 
taking into account that the eigenvalue occurs with mul- 
tiplicity 2s, we finally obtain the negativity 



M 



0,- 



s + 2(si-S2) 



2s + 1 



(11) 



The negativity is only determined by a single correlator, 
which is due to the fact that our state is SU(2)-invariant. 
At this stage, to obtain straightforwardly analytical 
results of negativity in the present ferrimagnetic system 
as well as to gain some essential physical insight into en- 
tanglement features like the intriguing high-temperature 
entanglement, we here focus only on a two-spin case first. 
In this case, the Heisenberg interaction can be written in 
terms of projectors as follows 



H = si 



S2-^[(S1 

2^ 



-S2f-S?- 
S2)'-|- 



s{s + 1) 



8 



1, 



s-1/2 



^P, 



+ 1/2- 



2 '""'" 2 
From Eq. (|12() . the eigenvalues are simply 



En 



s + 1 



El 



(12) 



(13) 
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FIG. 1: Negativity versus temperature for different s in the 
two-site system. 



The eigenvalue is just the expectation value of the corre- 
lator (si • S2) in the corresponding state. From Eqs. H13|) 
and (|ll|l . the negativity for the ground state and the first 
excited state are obtained as 



^fo 



1 



2s + 1 



, AAi = 0. 



(14) 



Clearly, the ground-state is entangled, while the first- 
excited state is not. In addition, we observe that the 
ground-state entanglement decreases as s increases, and 
finally vanishes in the limit of s — 00. Note that, the 
case of s = cxD corresponds actually to a classical spin, 
and thus, as expected intuitively, there exists no entan- 
glement between a quantum system and a classical one 
indeed. 

At finite temperatures, we need to know the correla- 
tor (si • 82} evaluated on the thermal state, which com- 
pletely determines the negativity. From Eq. I|13(l . we ob- 
tain readily the partition function and the correlator as 



Z = 2se^^ + 2(s + l)e~~ 



(15) 



and 



/ c \ 1 '5^ s(s + l)/i+M _££\ ._,„. 

(si-S2)=--^= ^{'' ' -' ')■ (16) 

Substituting Eqs. ifTB)) and ifT^ into Eq. flTT) . the nega- 
tivity is explicitly expressed as 

A/'=(^j^max{0,e^^^-2(s + l)e-^} (17) 

The numerical results of negativity versus temperature 
for different s are plotted in Fig. 1. As expected, for a 
fixed s, the negativity decreases monotonically as tem- 
perature increases, and vanishes when the temperature 
is equal or larger than a threshold value Tth. But more 




FIG. 2: Threshold temperature versus s. 




FIG. 3: Negativity versus temperature for different s in the 
four-site system. 



arrestingly, when s is large, the entanglement is rather 
robust against temperature. 

From Eq. I|17() , the threshold value of the thermal en- 
tanglement is found to be 



Tth = 



2s + 1 



21n(2s + 2)' 
from which we immediately have a striking result 



lim Tth = oo. 

S — *CXj 



(18) 



(19) 



The threshold temperature Tth can be arbitrarily high 
when s is large enough. A plot of Tth versus s is il- 
lustrated in Fig. 2, and we observe that the threshold 
temperature increases with the increase of s. 



It is now interesting to estimate some typical values 
in terms of Eq. (18) (with the unit J/fcs). The value 
of exchange constant J is of the order meV. For a tem- 
perature T ^ 127K, the corresponding s is estimated as 
s — 50, above which the thermal entanglement exists. In 
this case, the negativity is approximately 0.01 at temper- 
atures well below the above temperature. 

Finally, we wish to address briefly the negativity of 
the systems with a larger number of sites L. For the 
four-site case, we can still have exact results and plot the 
negativity versus temperature for different s in Fig. 3. 
It is seen that the negativity behaviors of the four-site 
case are qualitatively the same as those of the two-site 
case, namely, when s increases, the negativity at zero 
temperature decreases and the threshold temperature in- 
creases. As for even larger sizes L = 2N, due to the limi- 
tation of computational resource, we here can only make 
a qualitative analysis with the help of an approximate 
analytical result for the ground state energy obtained 
from the spin- wave theory JO. At zero temperature, 
the nearest spin correlator 1Z| (si ■ 82)0 = Eq{L)/2L = 
—s/2 — 5{s), where the subscript index denotes the 
ground state, and 5{s) > and approaches (1/4) in the 
limit s — > 00. Therefore, from Eq.(ll), the negativity 
JV = 25{s)/{2s -f 1) is non-zero at zero temperature. 
Also, for a small size chain, we note that there exists 
a energy gap (proportional to s) between the excited and 
ground states, so the negativity is mainly determined by 
the ground state contribution at finite temperatures well 
below the gap energy, and may be non-zero at high tem- 
perature for a very large s. 

In conclusion, by considering a simple ferrimagnetic 
chain model, we have found that the entanglement is 
rather robust against temperature. As the spin s in- 
creases, the threshold temperature for entanglement can 
be arbitrarily high for a small size chain. We hope that 
the present work motivates interests to investigate other 
physical systems which display high-temperature entan- 
glement. 
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